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Abstract 

We have calculated the density of states J\f(u>) of a dirty but homogeneous 

superconductor in a high magnetic field. We assume a dilute concentration 

of scalar impurities and find how N{oo) behaves as one crosses from the weak 

scattering to the strong scattering limit. At low energies, AT(lo) ~ to 2 for 

small values of the impurity concentration and scattering strength. When the 

disorder becomes stronger than some critical value, a finite density of states 

is created at the Fermi surface. These results are a consequence of the gapless 

nature of the quasiparticle excitation spectrum in a high magnetic field. 
PACS numbers: 72.20.My, 71.45.-d, 71.30.+h 
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I. MOTIVATION 



When an electronic system is placed in an external magnetic field, motion of electrons 
perpendicular to the field direction is confined to cyclotron orbits. The electronic energies 
are quantized in the form of a discrete set of Landau levels (LL's) separated by hu c , where 
u c = eH/mc is the cyclotron frequency. In high magnetic fields and at low temperatures, 
LL quantization leads to unusual behavior in many-body systems. Perhaps the most vivid 
example of such behavior is the fractional quantum Hall effect in a two-dimensional (2D) 
Fermi system subjected to a strong perpendicular magnetic field. There has also been 
much interest lately in the influence of LL quantization on properties of superconductors 
in a magnetic field. It was shown by Tesanovic et.al. [[jj], that standard Abrikosov-Gor'kov 
theory of type-II superconductors which neglects LL quantization, necessarily breaks down 
at high fields and low temperatures. The inclusion of LL's in the BCS description of the 
superconducting instability leads to reentrant behavior at high fields (much larger than 
the semi-classical upper critical field H C 2(T)) where the superconductivity is enhanced by a 
magnetic field [|IJ. As a consequence of the underlying LL structure, H C 2(T), or rather T C (H), 
develops oscillations near H C 2(0). Similar types of quantum oscillations, with the same 
origin, have been predicted in various other measurable quantities which are particularly 
pronounced in 2D systems 0. An important theoretical question in this context is the 
nature of the quasiparticle excitation spectrum at high magnetic fields. It is well known 
that the existence (or, typically, absence) of low-energy excitations in the superconducting 
state shapes the behavior of all superconducting thermodynamic and transport properties. 
The problem of quasiparticle excitations has been studied extensively in the low field limit 
(H > H c i(T)) in the context of an isolated vortex line. Scanning tunneling experiments 
have revealed the structure of bound states in the cores of isolated vortices at low fields, 
thereby confirming an early prediction of localized midgap states by Caroli, de Gennes and 
Matricon Only recently, was the problem of quasiparticle excitations in high magnetic 
field in the presence of a vortex lattice addressed by a number of groups At high 
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magnetic fields in extreme type-II superconductors vortices are closely packed (H < H c2 (T)) 
so that the isolated vortex picture from the low-field limit necessarily breaks down. In clean 
samples at low temperatures quasiparticles can propagate coherently over many unit cells 
of the vortex lattice. The solution of the BCS problem |J, where the coherent nature 
of quasiparticle excitations is fully accounted for, points to the qualitatively new nature 
of the energy spectrum at high fields: at fields near H C 2 this spectrum is gapless at a 
discrete set of points on the Fermi surface. This gapless behavior in 3D systems persists 
to a surprisingly low magnetic field ~ 0.5H C 2 ||, below which the gaps start opening up, 
and the system eventually reaches the regime of localized states in the cores of isolated 
vortices [[UJ. The gapless character of the excitation spectrum leads to algebraic behavior 
of various thermodynamic functions in the high field/low temperatures portion of the H- 
T phase diagram ||. In particular, it was shown by Dukan and Tesanovic (Tl[| that the 
recent observations of the dHvA-oscillations in the mixed state of A- 15 superconductors 
|i2| , [T3|| follow from the presence of a small portion of the Fermi surface containing gapless 
quasiparticle excitations, surrounded by regions where the BCS gap is large. This picture 
has received further support in the dHvA experiments on YNi 2 B 2 C | |14| |. 

In this paper we address the influence of the non-magnetic (scalar) impurities on the 
gapless superconducting state in high fields, by examining the behavior of the supercon- 
ducting density of states (DOS) in the presence of disorder. We are interested in a dirty 
but homogeneous superconductor for which the coherence length £ is much longer then the 
effective range £j mp of the impurity potential. Under this condition the order parameter in 
the mixed state is not substantially affected by the impurities and it is assumed to form a 
perfect triangular vortex lattice. In the opposite, inhomogeneous, limit the impurity poten- 
tial can pin the vortex making an otherwise perfect vortex lattice disordered. The problem 
of the superconducting instability and the quasiparticle excitation spectrum in the presence 
of randomly distributed vortex lines has been addressed so far only by Gedik and Tesanovic 
in the quantum limit, where electrons occupy only the lowest LL. It was demonstrated that 
the gaplessness of the excitation spectrum found in a pure superconductor is preserved even 
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when the positions of the vortex lines are completely random |15] . 

To understand the effect of disorder on the superconducting properties in the high-field 
limit near //^(O) we start with the normal state in this limit. This problem has been 
thoroughly investigated in the study of transport properties in high fields |[L6|| , where it was 
found that the electronic scattering with scalar point-like impurities randomly distributed 
through the sample leads to isotropic broadening by 5E = h/2r of the (LL's), where l/2r is 
the scattering rate due to the disorder. As long as the broadening SE -C hu c , the discreteness 
of the LL structure is preserved and the effect of the disorder on the superconducting state 
can be analyzed pertubatively starting from the results for the pure superconductor in high 
field. We present a Green's function pertubative approach to impurity effects in this regime 
based on the theory of superconducting alloys of Abrikosov and Gor'kov [[17]] This is a 
standard field-theory technique of treating the disorder in the superconducting state, and 
has been used extensively in the study of anisotropic superconductors such as heavy-fermion 
systems, p- or d-wave unconventional superconductors |21[ or in the study of superfluidity 
in 3 He films 0. 

This paper is organized as follows: In Sec. II we describe a simple model of electronic 
scattering from point-like scalar impurities randomly distributed in the superconducting 
material. In Sec. Ill we develop a self-consistent procedure for calculating the electronic 
self-energies in the framework of the Self-consistent Born (SCBA) approximation. We discuss 
the behavior of the superconducting density of states within this approximation. In Sec. 
IV we study the low-energy properties of the density of states in the self-consistent T- 
matrix approximation. In Sec. V we discuss the consequences of our results on various 
experimentally measurable properties of a superconductor in high magnetic field. 

II. DESCRIPTION OF THE MODEL 

We consider a 3-dimensional (3D) BCS-type, weak-coupling electronic system in a high 
magnetic field. We assume that the system in question is translationally invariant in the 
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normal state. The aspects associated with the true band structure of a particular material 
are not important for the issues that we are planning to address and are not discussed in 
this work. Nevertheless, all the results presented here can be modified to include band 
structure effects. For a weakly to moderately interacting system, one uses a simple short- 
range attractive BCS model interaction V(fi,f2) = — V5(ri,r^) between electrons. The 
dynamical origin of this interaction (i.e. whether mediated by phonons, charge density- 
fluctuations, etc.) is not important for our present purposes. We assume that the presence 
of impurities does not effect the effective electron-electron interaction so that it is still of the 
BCS form. Furthermore, the coupling constant V is taken to be independent of magnetic 
field. The mean-field (MF) Hamiltonian for this model system is: 

H = £ / *«(*0[^-HfcV + -A) 2 5 aP + Uapif) - g^ ■ H(r) - //]^(r)d 3 r + 

y"A(f)^{(f)^{(f)c/ 3 r + /i.c. (1) 

where ^ a (r) are the electron field operators for two spin components and fi is chemical 
potential. A is the vector potential due to the external field H, which is taken to be 
uniform everywhere in the system. Term — g/iecr ■ H describes the Zeeman splitting and 
U a p(r) = J2j U a p(r — R l ™' p ) is the random impurity contribution. A(f) = V < ^(r)W^(r) >, 
with < ... > denoting thermodynamic average, is the superconducting order parameter. 

For simplicity we consider a system in which Zeeman splitting is negligible, i.e. g ~ 0. 
The results are straightforwardly generalized to the g ~ 2 case [§|-|9|]. We assume that there 
are only non-magnetic (scalar) impurities present in the sample so that U a p(r) does not 
contain any spin-exchange terms. We consider dilute impurity concentrations for which the 
MF picture presented in this paper is valid. For a dirty but homogeneous superconductor 
in which the coherence length £ is much longer than the effective distance £j mp over which 
the impurity potential changes i.e., £/£j mp ^> 1, the superconducting order parameter A(r) 
in Eq. (1) is not affected by the impurities apart from its overall magnitude, and forms a 
perfect triangular Abrikosov lattice [18]. A generic example of such an impurity potential is 
an infinitely short-range 5-function potential of the form: 



UapW = E U{f-Ri) = £ U 5(f- Ri) (2) 

i i 

where Ri is the location of the zth impurity taken to be completely randomly distributed 
everywhere in the sample. The scalar scattering amplitude U a is assumed to be isotropic. 

It was shown in Ref. 5 that the unperturbed part of MF Hamiltonian (1) with U a p(r) = 
can be diagonalized in terms of the basis function of the Magnetic Sublattice Represen- 
tation (MSR) |TjJ, characterized by the quasi-momentum q perpendicular to the direc- 
tion of the magnetic field. The eigenfunctions of this representation in the Landau gauge 
A = H(—y, 0, 0) and belonging to the mth Landau level are: 



— exp (ik z () E ex P iy-^k 2 - ikq y b y ) 



—\ T T T 1 V"' U ^S I / j WV F V " q 

2 n 7l!v/?/ \ L x^y^z k 

exp [i(q x + —)x - l/2(y/l + qj + —l) 2 }H rn (^ + (q x + — )/). (3) 
where ( is the spatial coordinate and k z is the momentum along the field direction, a = (a, 0) 



and b = (b x ,b y ) are the unit vectors of the triangular vortex lattice, I = yhc/eH is the 
magnetic length and L x L y L z is the volume of the system. H m (x) is the Hermite polynomial 
of order m. Quasimomenta q are restricted to the first Magnetic Brillouin Zone (MBZ) 
defined by vectors Qi = {b y /l 2 , —b x /l 2 ) and Q 2 = (0,2a// 2 ). In the Landau gauge the 
Abrikosov order parameter can be written as: 

A(r) = A E ex P (in—n 2 ) exp (i2nnx/a — (y/l + nnl/a) 2 ) (4) 

n a 

where A is the overall BCS amplitude. The above form of the order parameter is taken to 
be entirely contained in the lowest LL of Cooper charge 2e. This is an excellent approxi- 
mation in the high-field regime QX|] . Normal and anomalous Green's functions for the clean 
superconductor in this representation can be constructed as: 

g{f,r'-u)= E 4> ntkz j(f^4>l >kzi? (P)G n (k z ,q;uj) 

n,k z ,q 

P(f, ?■ u) = £ K,-k z ,- q <^Kk*Ar')F*(k z , q; u) (5) 

n,k z ,q 
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where uo = kBT(2m + l)ir are the electron Matsubara frequencies. In writing (|5|) we have 
taken into account only diagonal (in Landau level index n) contributions to the Green's 
functions. (|5|). This is an excellent approximation in high magnetic fields where A/hu c <C 1 
and number of occupied Landau level n c is not too large. In this situation we can use the 
diagonal approximation JTJ, in which the BCS pairs are formed by the electrons belonging to 
the mutually degenerate Landau levels at the Fermi surface while the contribution from the 
Landau levels that are separated by hu c or more is included in the renormalization of the 
effective coupling constant (V — > V(H,T)) ||. In lower fields, where n c is a large number, 
the off-diagonal terms in (||) should be included on equal footing. It was shown in Ref. 
5 by numerically solving the BCS equations that the off-diagonal pairing does not change 
the qualitative behavior of the superconductor in a magnetic field as long as the magnetic 
field is larger than some critical field H critica i(T) (estimated to be ~ 0.5H C 2 at T pa in 
A- 15 superconductors): Placed in a high magnetic field and cooled to low temperatures the 
type-II superconductor has a gapless excitation spectrum. In the diagonal approximation 
the gapless branches in the spectrum can be found analytically as: 

E n (k z , q) = ±^/el(k z ) + \A nn (q)\* 

h 2 k 2 

e n (k z ) = + huj c (n + 1/2) - /i (6) 

where A nm (q) are the matrix elements of the Abrikosov order parameter ([|) in the MSR 
representation and can be written as: 

A (—1)'" b x 
A nm (q) = -7=— — r ==y2exp(iir—k 2 + 2ikq y b y - (q x + irk/a) 2 l 2 )H n+m [V2(q x + nk/a)l] 
Vz 2™ +m vn'm! & a 

(7) 

Once the off-diagonal pairing is included, the excitation spectrum cannot be written in the 
simple form ([]) and a closed analytic expression for the superconducting Green's function 
cannot be found. Nevertheless, the qualitative behavior of the quasiparticle excitations, 
characterized by the nodes in the MBZ, remains the same. The main role of the off-diagonal 
terms in Eqs. (5) for magnetic field strengths such that H > H C riticai is to renormalize 
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the slopes around the nodes. Once the magnetic field is lowered below H cri ti ca i } the gaps 
start opening up at the Fermi surface signaling the crossover to the low-field regime of 
quasiparticle states localized in the cores of widely separated vortices |L(J. In this paper 
we are interested in how the disorder affects the excitation spectrum in the gapless regime 
and, in particular, how the DOS behaves in this high-field (and low-temperature) gapless 
regime of a dirty superconductor. For this purpose, we argue that taking only the diagonal 
terms in Eqs. (5) will correctly capture the low-energy behavior of DOS for a wide range 
of magnetic fields (as long as H > H critica i) while considerably reducing the computational 
difficulties introduced by the off-diagonal terms. 

Proceeding in the frame described above, the "Fourier transformed" (in the quasimo- 
mentum space) Green's functions G^(k z , q; uj) and F®*(k z , q; uj) for the pure superconductor 
can be easily calculated as: 

-iuj - e n {k z ) 



G° n (k ZJ q;uj) 



u 2 + El(q,k z 



with E n (k z , q) and A nm (q) given by Eqs. (6) and (7), respectively. 



III. SELF-CONSISTENT BORN APPROXIMATION 

The equations of motion describing the mixed superconducting state in the presence of 
scalar impurities are: 

1 



iu + — [V r - - teA(r}} 2 + M - £ U(f- fy) 



Q(f, r'; cu) + A(f)J-'\f, r'; uj) = 5(f— r'), 

T ] {r, ?■ uj) - A*(f)£(f, ?■ uj) = (9) 



where normal and anomalous Green's functions incorporate the interactions between the 
electrons and impurities. We are interested in the quantities Q and averaged over the 
positions of the impurities. 
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In performing the average over the disorder we follow closely the diagrammatic technique 
developed by Abrikosov and Gor'kov in Ref. |I7| for the study of the superconducting alloys 
in zero field. First, we notice that for a dirty but homogeneous superconductor the quasi- 
momentum q is still a good quantum number, since the Abrikosov lattice is not affected 
by a short-ranged impurity potential (Q). Therefore, the Green's functions in (Q) can be 
expanded in terms of a complete set of eigenfunctions (H) as: 

g(r, r'- to) = £ E K <f)K (r')G(h, k 2 ; u) 
fci ki 

T\ry-u>) = j;E^(^W*i.fe«) ( 10 ) 

fel fe2 

where k = (q, k z ,n). 

The Hamiltonian for the interaction with impurities contains operator products 
Therefore, when an impurity is inserted into an electron line, two possibilities arise for each 
of the propagators Q, T and . These possibilities can be written in matrix form: 

G(x, x') -> G(x, y)a z G(y, x') (11) 

with the 2x2 matrix Q(x,x') defined in Nambu formalism, a z the Pauli matrix, and x = 
(r , r) , where r is the imaginary time. Taking into account the expansion (|10D this matrix 
equation can be rewritten in terms of its "Fourier" components as: 

G(h, k 2 ; u) = G°(h; u)5(h - k 2 ) + G°(h; 00) £ U klk3 (R t )G(h, k 2 ; 00) 

-F°(k Vl uj) £ U^ kl (Ri)F*(h,k2;uj), 

F*(k u k 3 ;uj) = F°*(kx;u)8(ki - k 2 ) + F°* (h; 00) £ U klk3 (R t )G(k 3 ,k 2 ;uj) 

k2,Ri 

+G°(-k 1 ;-uj) E U-. ka -k 1 (Ri)F*(fa,k 2 ;u) 

(12) 

with a similar equation for the function F(ki,k3',u). U kl k 2 {Ri) is the matrix element of 
the impurity potential U(f — B4) between two eigenstates (|3|). In high magnetic field we 
can assume that the scattering potential is weak compared to the separation between LL's 
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(given by hu c )- Under these circumstances electrons scatter into the states belonging to the 
same LL. Therefore, in solving Eqs. (12) we will neglect the inter-Landau level scattering 
that becomes important only at much lower fields. 

We first solve Eqs. ( |T2"D in the Born approximation by summing the diagrams that 
describe two consecutive electronic scatterings from the same impurity. Since the impurity 
atoms are randomly distributed through the system, we have to average expressions (|i~2"D 



over the position of each impurity. For a dilute concentration of scatterers with uncorrelated 
positions, we encounter two types of averages: 



< U klk2 (Ri)U k2ks (Ri) >$.— TliUl V r-i^ife^Li^bna^mnalft _ &) 

x y z v 



< U klk , 2 (Ri)U-k 3 -k 2 {.Ri) >R- "^0777 7^^ z i,fez3^i,g3^ni,n3'S'n 1 n 2 (9l ? , ?2) 

Li x LiyLi z y TIL 



(13) 



(14) 



where < ... >^ denotes the average over the impurity positions and n« is the impurity 
concentration. We find by inspection that (q± — q^) in Eq. ([13]) is a weak function 
of quasimomenta q\ — q 2 and LL indices rti and n 2 . For rti + n 2 > 5, to a very good 
approximation, S% (q) is independent of either q or n. The function S£ in2 (qi, (&) in Eq. 
( [TID is obtained as: 

^ ™«[«i.™2] (2k - 1)!! 

S nin2 {qi,q2) = V2 7oU\\] /("l-fc)(wl-fc)(gl)/faa-fc)(na-fc)(ga) ( 15 ) 

k=o y AK ) A 

where matrix elements f n n(q) = A nn (g)/A are calculated in (^) for the order parameter 
A(r) from the lowest LL of BCS pairs (|]). 

Using the above averages and the expressions (|]) we can bring the set of Eqs. ( jT2|) to 
the form: 

(iw - e n (k z ) - E N (tu)) G n (k z , q; u) + (A nn (q) + w)) F n *(fc„ «f; w) = 1 , 

(iu + e n (k z ) + E N (-u)) F^k z , q- u) + (A* nn (q) + E^(q, uj) G n (k z , q; u) = , 



(16) 
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The diagonal (normal) self-energy H N (u) and off-diagonal (anomalous) self-energy T,^ n (q;uj) 
can be expressed as: 

Z N (u) = n t U 2 ) £ G m (k z ,k;u) 

L x L y L z ^7Tl 

KM w) = rkPg . , ) n E k)F m (h, k; u>) (17) 

Li x Liylj z y TIL -, 

w m,k z ,k 



where the functions G m (k z , k; lo) and F m (k z , k; to) are found as a solution of Eqs. ([TBI) and 
can be written in the form: 

^ /, r \ -iuj-e m (k z ) 
& m {k z ,k;uj) 



^ + el l (k z ) + \A mm (kW 



F m (k z , k; u) = (i 8 ) 

^ + el(k z ) + \A mm (kW 



with 



iu = iu - £ (u), A mm (k) = A mm (k) + £ mm (/c; u) . (19) 

The sums over k z in Eqs. ( |TT| ) can be readily done so that Eqs. fliPP become self-consistent 
equations for the self-energies: 

| Tic r fX\ 1 f (jj 

u = u + To—— V — - — / dk . (20) 



and 



a ^ a r ° ^ V2(2fc-1)H , _ 
A nn (gJ = I\ nn {q) + 2^ (2k~)\\ Hn-k)(n~k){q> 

x ^f^)H)W^W (21) 

v=k A^k Fv J Ju 2 +\K uv {k)\ 2 



where iV(0) is the density of states at the Fermi level of the normal metal in zero field. 
T = TmJJlNiQ) is the scattering rate due to the disorder (defined in zero field). It is 
assumed that Y$/Ep <C 1 within the Born approximation. 

The exact self-consistent solution of Eqs. ( |2"0"| ) and ( f2T| ) is extremely difficult to obtain 



due to the coupling of matrix elements A nn (q) with A( n _k)(n-k)(<l) in Eq. (|2~I1) . Nevertheless, 
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this problem can be simplified if we notice that the behavior of A nn (q) for different LL indices 
is very similar around the gapless points in Eqs. (|6]) and differs considerably only in the 
regions in the MBZ that are gapped by large A. We are primarily interested in the low-energy 
behavior of the density of states and low-temperature thermodynamic properties of a dirty 
superconductor in a high magnetic field. These properties are governed by the quasiparticle 
excitations around nodes in A nn (q). Therefore, retaining only the k = term in Eq. (|2~Tf) 
represents a reasonable approximation in solving Eqs. ( P0f) and ([H]). Furthermore, we have 
demonstrated numerically that the sum of k ^ terms in (|2lD is < 10% of the k = term 
for small frequencies, while it is negligible for higher frequencies (of order A). Within this 
approximation A nn (q) = Af nn (q) so that Eqs. (^) and ( pT|) are combined as 



" ~ A + C iV(0) £ M Fn J " k ^2 + 1^)12 <22) 

where u = Co/ A and ( = r /A. The amplitude A = A(H,T,T) has to be determined from 
the self-consistent equation: 

A(f) =VTj2F(f,f;uj) (23) 

If we take the order parameter A(r) entirely in the lowest Landau level for the Cooper pair, 
Eq. (p3|) can be rewritten as 

LJ y L z 4r tin ^ + <k z )* + I A rm (q)\ 2 

where V(H, T) is the BCS pairing interaction that is, in general, renormalized by off-diagonal 
terms due to the coupling of electrons from LL's separated by hu c or more 0. One can 
think of V(H, T) as being chosen to reproduce the true self-consistent A(H, T) in a formalism 
that keeps only the diagonal terms. Its explicit form is easily computed in the A/hu c 1 
regime 0. Note that, once disorder is included, V(H,T) itself ha.s to be recomputed self- 
consistently. Here we ignore this complication on the grounds that a modest sacrifice in 
quantitative accuracy (for A/hu c < 1) is justified in the face of overwhelming numerical 
difficulty in determining self-consistent A(H,T) in presence of disorder. 
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The set of Eqs. and ( j2"4|) completely describes the effect of disorder on the supercon- 
ducting state within the Born approximation and enables us to calculate various physical 
quantities. The superconducting density of states in the presence of impurities is defined as 

lr 1 

M a {uj) = — Sm / dfG(r,f;iu) = — — — Sm G n{k z ,q;iu)\ iuj=ul+i5 (25) 

7T J 7T IJt. 1j<ii±j 7 , 

* y * n,k z ,q 



where G(k z , <f; ico)\iu=w+is is given by expression (|18D in which the analytic continuation to 
real frequencies is performed. Equation (p2|) is an implicit equation from which u = u[u/A] 



is to be calculated. Once u is known, the density of states (25) can be obtained as 



In Fig. 1 we plot Af 8 (u)/N(0) for several values of the parameter ( = Tq/A when n c = 10. 
Two kinds of behavior are present: For ( < 0.9 the superconducting density of states vanishes 
at the Fermi level and J\f a (uS) ~ uo 2 for small u. This is the same behavior as one finds in a 
pure system, just that the coefficient in front of uj 2 is increased from its clean system value. 
When C, > 0.9, a finite density of states is created at the Fermi level, although of course it is 
smaller than in the normal state. In this regime J\f s {w) ~ A/" s (0) + const, x uj for small u. At 
higher energies, we observe that the peak in the density of states located at to /A fa l/v2 
in the clean system is reduced and broadened by disorder. As the impurity concentration 
(measured by the parameter Q increases, the peak eventually disapeares. Note, though, 
that our calculation might be less accurate at the higher energies (of the order A) due to the 
number of approximations that are, as explained above, strictly applicable only at the low 
energies. Furthermore, the true behavior of the peak can be investigated only if A = A(r, T) 
is found from the self-consistent equation (p4|). 

IV. T-MATRIX APPROXIMATION 

The self-energies S,^ n (g;u;) and S^ n (g;u;) of the superconducting system obeying Eqs. 
(PS|) are closely related to the diagonal (with respect to the Magnetic Translation Group 
basis) T-matrix elements in a single-site approximation as 
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T$ n (k;u>)=n i <T u (k,lr,u>) > n% 
^ n (k;u) = -n l <T 12 (k,k;u J )> Ai (27) 

where T* 3 '(fci, A^;^) are the coefficients in the T-matrix expansion over the complete set of 
eigenstates k = (k z , k, n) 

T n (f,f » = £ fcl (rV£ 2 (r ')T n {hMu) 

ki,k2 

T 12 (f,f » = J2 (t>-k 1 (f f )K(r , )T n (k 1 ,k 2 ;uj). (28) 
The 2x2 T-matrix T(r,r';uj) obeys the Lippmann-Schwinger equations 

T(f,f';uj) = U(r}S(f—f')a z + / dfiU(f)a z Q(f, ri;o;)T(fi,f ';o;) (29) 



where matrix elements are given by (|T0|) and Z7(r) is the impurity potential (g). As 
in SCBA, we neglect the inter-Landau level scattering on the basis that in high magnetic 
field the scattering potential is much weaker than hu c . After averaging over the impurity 
position, Eqs. ( p9| ) for the dirty but homogeneous superconductor reduce to 

T nl( k z, k; uj) = U + Uq l * — G m(Qz, q; uj)T^ n (q z , q; k z , k; u) 

x V 2 V q z ,q,m 

+U by r S nm(k,q)F m (q z ,q;u)T^ n (q z ,q-,k z ,t,^) (30) 

x y * v q z> q,m 

and 

Tnn(h, k; u) = -Uo L J^J—T^j Yl S nm(k, q)F^(q z , q; u)T^ n (q z , q; k z , k; u) 

x V z v q z ,q,m 

+U T T by T -j=i ^ G m(-(lz,-q;-^)T^ n (q z ,q;k z ,k;^)) (31) 

y *• v q z ,q,in 

where G m (q z , q; to) and F m (q z ,q;uj) are written in the form (|I8"D and S^ m (k,q) is given by 
formula (|1|). Around the gapless points of the excitation spectrum ([]) the second term in 
equation (|30| ) is very small compared to the first term. This can be deduced if one inspects 
expressions flTSP and fl3"T| ) around the nodes of the excitation spectrum ([]) in the MBZ. 
Keeping this conclusion in mind, equations fl30"l) and (|3lf) can be solved as 
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and 



rpllt \ L x L y L^l^qz,q,rn G m(<lz,q,m) ^ 

~~ — 2 — f — & ~ — ^ 



^^lfnn(k) J2q z ,q,m frnrn(q)F m (q z , q\ id) 



T*l(k; u) = _ Ww y- — - — __ (33) 

where we used the explicit form ( |i5|) for S^ m (k, q) (taking only the k = term in Eq. (|15"D, 
see the discussion in previous section). With the help of definitions ([19]) and (|2T|), the set 
of equations (|32| ) and (|33]) can be brought into the form 



M = X + ^ r =V2 ( 34 ) 



<? + [En ^Cw Jdqu/y/u* + l/^l 2 ] 

where £ = T/A and w = u)/A. 

Disorder is characterized with two parameters: T = nj/iV(0)7r = (ni/n)Ep, which 
measures the concentration of impurities rij relative to the electron density n, and c = 
l/irN(0)Uo, which measures the strength of the scattering potential. The normal state in- 
verse scattering rate l/2r is found by taking f n n(o) = in (|34] ) and letting to — > 0. This 
procedure yields l/2r = T/(l + c 2 ), the result first obtained in Ref. 17 in the study of the 
transport properties of normal metals in a high magnetic field. The weak scattering limit is 
approached when c 2 is much larger than the second term in the denominator of expression 
([34]), while the strong scattering limit is achieved when c 2 = 0. In the strong scattering 
limit, the approximation in which the inter-Landau level scattering is neglected eventually 
becomes unphysical, unless the magnetic field is so high that only the lowest LL is occupied. 

As in previous section, the superconducting density of states in the presence of impurities 
is found from Eq. (p6| ) once the solution u = u[u/A] of Eq. (|34]) is found. Figures 2, 3 
and 4 show how the superconducting density of states J\f s (uj)/N(0) behaves as a function of 
the energy parameter u/A as one crosses from a weak scattering limit, c = 1.0, (Fig. 2) to 
a strong scattering limit, c = 0.0, (Fig. 4). For each value of c (measuring the scattering 
strength) we present how the density of states changes as the impurity concentration, i.e., 
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parameter ( increases. There are two types of behavior present in figures. When ( < c 2 the 
density of states vanishes at the Fermi level with J\f s (u) ~ uj 2 for small u, the coefficient 
in front of uj 2 being increased from the clean system value. When ( ^ c 2 , the density of 
states still vanishes at the Fermi level, but jV s (uj) ~ uj for small uj. Further increase of the 
concentration rij such that ( > c 2 creates a finite density of states jV s (0) at the Fermi level. 
In the strong scattering limit, c< 1, the superconducting density of states is finite at the 
Fermi level for any non-zero concentration of impurities. In this limit, j\f s (0)/N(0) ~ 2(7/A) 



where 7 = Ay£/2 for £ <C 1- Furthermore, below uj/A = y/£ in the strong scattering limit, 

a peak is observed centered at zero energy. This peak suggests formation of a quasi-bound, 
resonant state which is analog of a Shiba state formed in the energy gap of a conventional 
s-wave superconductor, as a result of multiple scattering off a magnetic impurity [pOR . As 
one moves away from the strong scattering limit, this zero-energy peak disappears. At 
higher energies, we find a similar behavior of the density of states to the one observed in 
the Born approximation of the previous section: The peak, located at uj/A « 1/V2 in 
the clean system, is reduced and broadened as the impurity concentration (measured by 
parameter () increases. Also, this peak is slightly shifted to higher uj/A as ( increases, 
suggesting a stronger reduction in BCS amplitude A = A(T, T) than what is found in the 
Born approximation. 



In this paper we have analyzed the influence of a dilute static disorder on superconducting 
properties in a high magnetic field. We considered a dirty but homogeneous superconductor 
for which the order parameter A(r) is not influenced by the presence of the impurities and 
still forms the perfect Abrikosov triangular lattice. We considered the weak-scattering limit 
within a self-consistent Born approximation while the strong-scattering limit was treated 
within a T-matrix approximation for superconducting self-energies. 

We found that for small impurity concentrations and weak scattering potentials the 




V. CONCLUSIONS 
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superconducting density of states behaves as E 2 for small energies E, the same behavior 
as that found for the pure superconductor in a magnetic field ||. When disorder becomes 
stronger than some critical value, a finite density of states (but still smaller then the normal 
state value) is created at the Fermi level. The finite superconducting density of states at the 
Fermi level signals the broadening of gapless points into gapless regions in the MBZ. It is 



interesting to mention that this behavior is similar to that of dirty superfluid He films p2 



and somewhat similar to the behavior of the density of states in anisotropic heavy-fermion 
superconductors plfl . 

The experimental property of a superconductor in which the absence of a quasiparticle 
gap over some region of the Fermi surface will be most obviously felt is the specific heat. In 
a clean system in high magnetic field at low temperatures c v ~ AT 3 , where A is the field 
dependent coefficient ||. In a dirty but homogeneous superconductor, instead of the T 3 
law one finds linear behavior at low temperatures with the coefficient reduced by the factor 
~ 2 (7/ A) from the normal state value. Detailed measurements of heat capacity at very 
low temperatures and high magnetic fields are not yet found in the literature. We propose 
a class of A- 15 superconductors as good candidates in which the linear temperature law 
of heat capacity at high magnetic field can be discovered. These systems have experimen- 
tally accessible upper critical fields and are clear examples of materials for which the LL 
quantization in high fields plays an important role |pT|-[T3|| . 

This work has been supported in part by the NSF Grant No. DMR-94 15549. 
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FIGURES 

FIG. 1. Quasiparticle density of states Af 8 (u) / 'N (0) vs reduced energy u>/A in the Born ap- 
proximation as a function of parameter £ = To/ A. 

FIG. 2. Quasiparticle density of states M s (u)/N(0) vs reduced energy to/ A when c = 1.0, as a 
function of the impurity concentration parameter £ = F/A. 

FIG. 3. Quasiparticle density of states Af s (u)/N(0) vs reduced energy uj/A when c = 0.5, as a 
function of the impurity concentration parameter £ = T/A. 

FIG. 4. Quasiparticle density of states Af s (u)/N(0) vs reduced energy lj/A when c = 0.0, as a 
function of the impurity concentration parameter ( = T/A. 
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